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ABSTRACT 

We use i?-band CCD linear polarimetry collected for about 12 000 background field stars in 46 fields 
of view toward the Pipe nebula to investigate the properties of the polarization across this dark cloud. 
Based on archival 2MASS data we estimate that the surveyed areas present total visual extinctions 
in the range 0™6 < Ay < 4™6. While the observed polarizations show a well ordered large scale 
pattern, with polarization vectors almost perpendicularly aligned to the cloud's long axis, at core 
scales one see details that are characteristics of each core. Although many observed stars present 
degree of polarization which are unusual for the common interstellar medium, our analysis suggests 
that the dust grains constituting the diffuse parts of the Pipe nebula seem to have the same properties 
as the normal Galactic interstellar medium. Estimates of the second-order structure function of the 
polarization angles suggest that most of the Pipe nebula is magnetically dominated and that turbulence 
is sub-Alvenic. The Pipe nebula is certainly an interesting region where to investigate the processes 
prevailing during the initial phases of low mass stellar formation. 

Subject headings: ISM: clouds — ISM: individual objects: Pipe nebula — Stars: formation — Tech- 
niques: polarimetry 



1. INTRODUCTION 

The relatively low Galactic star formation efficiency 
(SFE, defined as the fraction of a molecular gas mass 
that is converted into stars) is one fundamental con- 
straint on the global properties of star formation. In 
our Galaxy the SFE is observationally estimated to be 
of the order of a few percent when whole giant molecu- 
lar complexes are considered. For instance, the detailed 
stud y of the Taurus molec ular cloud complex conducted 
by IGoldsmith et al.M2008| ) provided a SFE between 0.3 
and 1.2%. Magnetic fields and supersonic turbulence are 
two mechanisms that are commonly invoked for regu- 
lation of such small SFE. Magnetic fields may regulate 
cloud fragmentation by several physical processes (e.g., 
moderating the infalling motions on the density peaks, 
controlling angular momentum evolution through mag- 
netic breaking, launching jets from the near-protostellar 
environment, etc). On the other hand, it is known that 
turbulence may play a dual role, both creating overdensi- 
ties to initiate gravitational contraction or collapse, and 
countering the effects of gravity in these overdense re- 
gions. The respective rules of magnetic fields and inter- 
stellar turbulence in regulating the core/star formation 
process are, however, highly controversial. For instance, 
some authors opine that magnetic fields are a bsolutely 
dominant in the s t ar formation pro cess (e.g., |Tassis fc| 
Mouschovias||2005| jGalli et al. 2006), while others sup- 
port that super- Aifvenic turbulence provides a good de- 
scription of molecular cloud dynamics, and that the aver- 
age magnetic field strength in those clouds may be much 
smaller than required to support them against the grav- 
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itational collapse (see Pad oan et al.||2004 and references 
therein). 

The Pipe n ebula, a massive fila mentary cloud com- 
plex (10 4 M0, |Lombardi et al.|2006|) l ocated at the solar 
vicinity (145 pc, Alves & Franco 2007) which presents an 
apparently quiescent nature, seems to be an interesting 
place to look for some answers on the physical processes 
involved in the collapse of dense cloud cores and how 
they evolve until stars are formed. 

Optical images of the Pipe nebula (see for instance the 
wonderful high quality image obtained by Stephane Guis- 
ard for the Giga Galaxy projectF]) or the dust extinction 
map obtained by Lombardi et al 



(2006), show that this 



map 

complex compris es many dark core s and sinuous dark 
lanes. Although |Alves et al] ( |2007[ ) have identified 159 
dense cores with estimated masses ranging from 0.5 to 
28 M© all over the entire Pipe nebula, the only known 
star forming active site in this nebula seems to be re- 
stricted to its northwestern extreme (in galactic coordi- 
nates), the densest part of the complex associated with 
the dark cloud Barnard 59 (B59), which corresponds to 
only a small fraction of the entire cloud mass. Actually, 
an embedded cluster of young stellar objects within B 59 
has been revealed by infr ared images obtaine d with the 
Spitzer Space Telescope (Brooke et al. 2007). The ap- 
parently low efficiency in forming stars o bserved for this 
cloud comp lex (only ~0.06% according to |Forbrich et al. 
2009, |2010| , suggests that the Pipe nebula is an example 
ol a molecular cloud in a very early s tage of star forma- 



tion. Indeed, in our previous paper (Alves et al. 2008 
hereafter Paper I) it has been suggested that the Pipe 
nebula may present three distinct evolutionary stages, 



1 http://www.gigagalaxyzoom.org 
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Fig. 1. — I dentification of the observed 46 lines-of-sight overplotted on the dust extinction map of the Pipe nebula obtained by Lombardi 
|et al.|p006| . The small squares roughly indicates the observed CCD field of view, which in our case corresponds to about 12' x 12'. The 
large retangles demarcate the areas detailed separately in Figs.[5] to [9] (colored version of this and of the other figures are available in the 
online version of this paper). 



being the B 59 region the most evolved of them while 
the opposite extreme of the cloud (the bowl) would be in 
the earlier stage. This suggestion seems to be reinforced 
by the recent Spitzer census of sta r formation activity 
performed by Forbrich et al. (2009), who detected only 
six candidate young stellar objects (YSOs) outside the 
B 59 region, four of them located in the "stem" of the 
Pipe, none having been detected in the bowl. Moreover, 
the youthfulness of the YSOs in B 59 is corro borated by 
the results obtained by Covey et al. (2010), who esti- 
mated a median age of about 2.6Myrs to the candidate 
YSOs found in B 59. Interestingly, they suggest that this 
population may be older than the well studied ones in 
Cha maeleon, Taurus, and p Ophiuchus, respectively. 

In |Paper~T| we described the global polarimetric prop- 
erties of the Pipe nebula as obtained from mean values of 
polarization degree and dispersion in polarization angles 
calculated for stars having P/o~p > 10. In the present pa- 
per we introduce the details of our data sample collected 
for 46 CCD fields, which are exactly the same as the one 
used in the previous work, and analyse the polarimetric 
properties of the Pipe nebula at core scales. In order to 
increase the statistical sample for each investigated field, 
we were less strict in our selection criteria accepting stars 
with P/op > 5. 



2. OBSERVATIONS 

2.1. Data acquisition and reductions 

The polarimetric data were collected with the 1.6 m 
and the IAG 60 cm telescopes at Observatorio do Pico 
dos Dias (LNA/MCT, Brazil) in missions conducted from 
2005 to 2007. These data were obtained with the use of 
a specially adapted CCD camera to allow polarimetric 
measurement s — for a suitable d escription of the po- 
larimeter see Magalha es et al.| ( |1996| ) . i?-band linear po- 
larimetry was obtained tor 46 fields (with field of view of 
about 12' x 12' each) distributed over more than 7° (17 pc 
in projection) covering the main body of the Pipe nebula. 
The observing lines-of-sight were visually selected from 
inspection of the IRAS 100 /im emissi on image of the 



Pipe nebula prior to the publication by |Lombardi et al 
( 2006 ) of the dust extinction map of this cloud complex. 



In our selection we chose directions toward high dust 
emission as well as some directions pointing to positions 
presenting lower emission but close to the main body of 
the complex as de fined by the 100 /im image. After that, 
Alves et al.| ( |2007| published their list of dense cores and 
some of our selected fields turned out either to completely 
include one of these cores or part of its outskirts. In Fig.[l] 
the observed lines-of-sight are overplotted on th e dust ex- 
tinction map of the Pipe nebula obtained by |Lombardi| 
et al. (2006). The small squares roughly indicates the 
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TABLE 1 

Observed zero polarization standard stars. 







ITurnshek et al.l 


ISchmidt et al.l 


this work 


HD 


V 


V v (%) 


V v (%) 


(%) 


12021 


8.85 




0.078 (0.018) 


0.106(0.037) 


98161 


6.27 


0.017(0.006) 




0.028(0.041) 


154892 


8.00 


0.050(0.030) 




0.027(0.041) 


176425 


6.23 


0.020 (0.009) 




0.031 (0.017) 


BD+28 4211 


10.51 




0.054 (0.027) 


0.066 (0.025) 



TABLE 2 

Observed high polarization standard stars. 



HD 



V 



ITurnshek et al.l 

W (%) 0"H 



IHeilesI 

p (kj—^o (°) 



this work 

p r (%) e (°) 



110984 
111579 
126593 
155197 
161306 
168625 
170938 
172252 



8.95 
9.50 
8.50 
9.20 
8.30 
8.40 
7.90 
9.50 



5.70(0.01) 
6.46(0.01) 
5.02(0.01) 
4.38(0.03) 



91.6 
103.1 

75.2 
103.2 



5.19 (0.11) 
6.21 (0.17) 
4.27(0.10) 
3.99 (0.08) 
3.69 (0.09) 
4.42 (0.20) 
3.69 (0.20) 
4.65 (0.20) 



90.6 
103.0 

77.0 
103.9 

67.5 

14.0 
119.0 
148.0 



5.21(0.19) 
6.11(0.09) 
4.65 (0.11) 
3.98 (0.07) 
3.59 (0.23) 
4.23 (0.07) 
3.62(0.11) 
4.38 (0.14) 



91.4 
103.1 

74.9 
105.2 

67.9 

14.9 
118.8 
147.7 



areas covered by the observed frames. 

When in linear polarization mode, the polarimeter in- 
corporates a rotatable, achromatic half- wave retarder fol- 
lowed by a calcite Savart plate. The half-wave retarder 
can be rotated in steps of 22? 5, and one polarization 
modulation cycle is covered for every 90° rotation of 
this waveplate. This arrangement provides two images of 
each object on the CCD with perpendicular polarizations 
(the ordinary, / , and the extraordinary, / e , beams). Ro- 
tating the half-wave plate by 45° yields in a rotation of 
the polarization direction of 90°. Thus, at the CCD area 
where f Q was first detected, now f e is imaged and vice 
versa. Combining all four intensities reduces flatfield ir- 
regularities. In addition, the simultaneous imaging of 
the two beams allows observing under non-photometric 
conditions and, at the same time, the sky polarization is 
practically canceled. Eight CCD images were taken for 
each field with the polarizer rotated through 2 modula- 
tion cycles of 0°, 22? 5, 45°, and 67? 5 in rotation angle. 

Among the 46 sky positions, twelve were observed at 
the IAG 60 cm telescope. At this telescope the integra- 
tion time was set to 120 seconds and 5 frames were col- 
lected and co-added for each position of the half-wave 
plate (totalizing 600 seconds per wave plate position). 
The remaining 34 fields were observed at the 1.6 m tele- 
scope, where the integration time for most of the ob- 
served positions was also set to 120 seconds, being that 
only one frame was acquired for each position of the half- 
wave plate. In order to have almost the same field of 
view, the latter telescope was provided with a focal re- 
ducer. 

The CCD images were corrected for read-out bias, zero 
level bias and relative detector pixel response. After 
these normal steps of CCD reductions, we identified the 
corresponding pairs of stars and performed photometry 
on them in each of the eight frames of a given field using 
the IRAF DAOPHOT package. From the obtained file 
containing magnitude data, we calculate the polarization 
by use of a set of specially developed IRAF tasks (PCCD- 



PACK package; |Pereyra||2000| ). This set includes a spe- 
cial purpose FORTRAJN routine that reads the data files 
and calculates the normalized linear polarization from 
a least-square solution, which yields the degree of linear 
polarization (P), the polarization position angle (6 , mea- 
sured from north to east) and the Stokes parameters Q 
and U, as well as the theoretical (i.e. the photon noise) 
and measured errors. The latter are obtained from the 
residuals of the observations at each waveplate position 
angle (ipi) with respect to the expected cos 4^ curve. 

Zero polarization standard stars were observed every 
run to check for any possible instrumental polarization, 
which proved to be small as can be verified by inspection 
of Table[l] The reference direction of the polarize r was 
determined by ob serving polarized standard stars (Turn- 
|shek et al.| [1990 ) , complemen ted with po larized stars 
from the catalogue compiled by Heiles ( 2000 ) . For all ob- 
serving seasons, the instrumental position angles showed 
a perfect correlation with the standard values (see Ta- 
ble^, and the expected uncertainty of the zero point for 
the reference direction must be smaller than 1-2°. 



2.2. Results 

Our final sample contains 11948 stars, being that 
9 777 of them have P/op > 5, where <jp means the 
largest between the theoretical and measured errors, that 
is, about 3 200 stars more tha n the ones used in the 
analysis conducted in |Paper l| which limited the sam- 
ple to stars presenting P/ap > 10. A search in the 
archival Two- Micron All- Sky Survey (2MASS), avail- 
able on-line http : / / irsa . ipac . caltec . edu| identified 
11588 objects that could be associated to our observed 
stars, and the JHK S photometry was retrieved for them. 
For the remaining we usually failed to associate a 2MASS 
object either because none was found inside the searched 
box (up to about 8 times the typical rms error of our as- 
trometric solution) or, in case of the most crowded fields, 
because the same 2MASS object could be assigned to 
more then one of ours. 
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Fig. 2. — Distribution of the polarimetric errors as a function 
of the J magnitude as retrieved from the 2MASS catalogue. The 
distribution shows characteristics of estimated errors dominated by 
photon shot noise. 

Figure[2] gives the distribution of the estimated polari- 
metric errors, op, as a function of the J2MASS magni- 
tude, for the observed stars. This figure shows that most 
of our stellar sample has magnitude within the interval 
10 m < J2MASS < 15 m and polarimetric error given by 
dp < 0.5 %. The obtained distribution suggests that the 
uncertainties are dominated by photon shot noise, as ex- 
pected for a sample collected with fixed exposure time. 

The distribution of obtained degree of polarization for 
stars with P/op > 5 (Fig. [3] - bottom panel) shows a 
surprising result: several stars present degree of polar- 
ization larger than 15%. Indeed, 6 objects present de- 
gree of polarization slightly larger than 19%. As far as 
we know, these are the largest polarization produced by 
dichroic extinction ever obser ved — in the stellar polar- 
ization catalogue compiled by Heiles ( 2000 ) one find only 
four stars, out of more than 9 000, with degree of polar- 
ization higher than 10%, being that the highest of them 
equals to 12.47%. 

The distribution of obtained polarization angles for 
stars with P/op > 5 (Fig. [3]- top panel) shows a large 
concentration of values around 6 « 180° (in equatorial 
coordinates), that is, 76% of this subsample has polar- 
ization angles between 160° and 10°, clearly indicating 
a high large scale homogeneity in dust grain alignment 
(and supposedly in the geometry of the magnetic field) 
all over the whole region. Since the main axis of the 
Pipe's stem is almost in line with the west-east direc- 
tion, it clearly indicates that the observed polarization 
vectors are mainly pe rpendicu larly aligned to the longer 
axis of the cloud (see Paper I). 

It is interesting to compare the distributions for stars 
having P/op > 5 with the ones having P/op > 10, also 
shown in Fig.[3J and used in our previous work. Visu- 
ally one can attest that the distribution of polarization 
angles for both samples are quite similar. Indeed, the 
mean value of both distributions differ by only 2° and 
the standard deviation of the sample having P/op > 5 
is about 2? 5 wider than the one for P/op > 10. On the 
other hand, as expected for a sample with uncertainties 
dominated by photon shot noise, most of the stars in- 
cluded by the less tight signal-to-noise condition are the 
ones presenting smaller degree of polarization. 

The high values of polarization obtained in our survey 
is the result of differential extinction produced by inter- 
stellar dust grains on the incoming/background stellar 
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Fig. 3. — The obtained distributions for the 9 777 stars with 
P/o~p > 5 and for the 6 582 stars with P/ap > 10, darker (salmon) 
and lighter (yellow) histograms respectively. Top panel: distribu- 
tion of the observed polarization angles; Botton panel: distribution 
of the estimated degree of polarization. 

radiation. It is assumed that a large fraction of those 
grains are aligned. The nature of such alignment is still 
a matter of debate (see |Lazarian||2003| for a comprehen- 
sive review on this subject), however, it is widely believed 
that the dominant process responsible for the alignment 
involves interaction between the spin of the dust particles 
and the ambient magnetic field, as originally proposed by 
Davis & Greensteinl (|1951|). 



3. DATA ANALYSIS 

3.1. Mean Polarization 

Apart from the difference in the a dopted signal-to- 
noise, the results described in Pap er I| and the ones pre- 
sented in this work were obtained in the way described 
below and introduced in Table[3] which gives, for each 
observed field, coordinates (right ascension and declina- 
tion), associ ated dark core, whe n available, from the list 
compiled by |Alves et al. (2007), the number of stars for 
which we estimated polarization, and the number of stars 
with P/crp > 5, in columns 2 and 3, 4, 5, and 6 respec- 
tively. 

Mean polarization and polarization degree were esti- 
mated for each of the observed area adopting a procedure 
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TABLE 3 

Mean .R-band linear polarization and extinction data for the 46 observed fields in the Pipe nebula (see text for 

explanation on the columns). 



Field 


alpha (J2000) delta 
h m s o / // 


Core a 


observed 
stars 


stars with 
P/a > 5 


n b 
Vmean 

O 


O 


stars 


(P) 

% 


SP 

% 


Q{p) 

o 


AO 

o 


A v 
mag 


5A V 
mag 


01 


1710 28 


-27 22 49 


06 


273 


174 


9.7 


10.65 


168 


2.45 


0.82 


3.1 


9.65 


2.54 


1.25 


02 


171152 


-27 03 49 


— 


165 


64 


12.9 


15.04 


63 


1.82 


0.85 


9.8 


13.72 


0.63 


0.35 


03 


171121 


-27 24 46 


12 


62 


23 


28.9 


24.93 


21 


1.78 


0.96 


48.8 


23.75 


3.88 


2.00 


04 


1710 55 


-2744 26 


— 


400 


211 


15.9 


13.08 


202 


1.63 


0.79 


14.8 


12.38 


1.59 


0.74 


05 


1712 30 


-27 20 42 


14 


137 


50 


0.9 


16.82 


47 


2.00 


1.07 


179.7 


15.67 


2.53 


1.09 


06 


1712 01 


-27 37 06 


08 


271 


91 


128.7 


40.79 


87 


0.62 


1.36 


149.8 


40.21 


2.33 


1.19 


07 


1713 53 


-2712 33 


— 


206 


114 


167.3 


12.75 


108 


1.27 


0.61 


169.3 


11.74 


2.06 


0.97 


08 


1713 34 


-2745 46 


— 


807 


486 


7.6 


8.80 


464 


1.76 


0.63 


9.2 


7.40 


1.58 


0.78 


09 


1714 52 


-27 20 55 


21 


189 


132 


6.3 


10.76 


54 


3.44 


0.98 


3.5 


10.04 


3.03 


1.49 


10 


1715 25 


-2718 09 


— 


191 


174 


4.2 


7.72 


164 


3.43 


1.23 


4.4 


7.29 


3.33 


1.67 


11 


171515 


-27 33 38 


20 


135 


130 


170.8 


6.76 


127 


3.17 


1.14 


171.3 


6.21 


2.32 


1.17 


12 


1716 20 


-27 09 32 


25 


198 


189 


176.5 


4.87 


179 


4.46 


1.01 


176.8 


4.32 


2.72 


1.22 


13 


171712 


-27 03 06 


27 


199 


196 


0.9 


5.76 


190 


4.03 


0.69 


1.3 


5.62 


2.03 


0.94 


14 


1716 05 


-273138 


23 


282 


280 


176.0 


4.40 


267 


4.34 


0.93 


176.1 


4.11 


2.82 


1.42 


15 


1718 27 


-26 47 50 


31 


272 


254 


7.7 


10.90 


243 


2.18 


0.78 


5.3 


10.63 


2.32 


1.16 


16 


1718 48 


-271136 


— 


382 


319 


2.2 


4.92 


301 


2.39 


0.50 


2.6 


3.52 


1.98 


0.98 


17 


1719 36 


-26 55 23 


33 


214 


210 


171.6 


5.91 


201 


2.58 


0.76 


171.8 


5.58 


2.07 


0.91 


18 


17 20 49 


-26 53 08 


34/40 


327 


201 


164.2 


9.59 


198 


4.63 


1.68 


163.0 


8.60 


3.64 


1.81 


19 


17 22 43 


-26 39 25 


— 


368 


323 


167.7 


7.06 


306 


2.31 


0.79 


167.2 


6.35 


1.86 


0.94 


20 


17 24 03 


-26 20 35 


— 


451 


400 


32.2 


9.86 


383 


2.05 


0.73 


32.5 


9.23 


1.49 


0.68 


21 


17 2148 


-271815 


— 


260 


218 


179.1 


8.73 


209 


1.96 


0.38 


0.8 


8.05 


2.33 


1.16 


22 


17 22 38 


-27 04 14 


41/42 


102 


97 


158.1 


5.35 


92 


4.39 


1.33 


157.0 


4.82 


3.03 


1.50 


23 


17 2713 


-25 07 27 


— 


513 


412 


7.1 


9.44 


391 


2.37 


0.74 


7.1 


8.79 


2.03 


0.99 


24 


17 26 25 


-25 58 09 


— 


748 


520 


175.1 


7.61 


506 


2.59 


0.76 


174.2 


6.52 


2.12 


1.05 


25 


17 28 07 


-25 29 52 


— 


511 


461 


174.2 


5.58 


437 


3.27 


0.74 


173.9 


4.91 


2.17 


1.08 


26 


17 25 40 


-26 43 09 


48 


247 


126 


142.3 


32.79 


120 


1.99 


1.49 


152.9 


32.52 


2.21 


1.26 


27 


17 25 28 


-27 03 29 


— 


197 


185 


143.2 


11.43 


183 


3.39 


1.30 


142.2 


11.19 


2.21 


1.10 


28 


17 3018 


-25 09 50 


— 


254 


240 


172.9 


5.86 


231 


5.61 


1.16 


173.3 


5.35 


3.18 


1.59 


29 


17 29 14 


-25 55 44 


-70 


98 


95 


160.2 


5.90 


91 


5.43 


1.52 


160.8 


5.46 


3.59 


1.76 


30 


17 2812 


-26 2110 


56 


94 


93 


160.6 


2.54 


90 


6.64 


2.17 


160.6 


2.42 


2.79 


1.39 


31 


172712 


-26 42 59 


51 


284 


271 


155.1 


5.12 


258 


4.99 


1.59 


156.1 


4.61 


3.47 


1.73 


32 


17 27 24 


-26 56 50 


47 


143 


139 


164.1 


4.58 


135 


6.20 


1.54 


164.5 


3.92 


3.81 


1.88 


33 


17 32 09 


-25 2418 


91 


329 


313 


169.9 


4.04 


301 


8.10 


1.30 


169.3 


3.22 


4.38 


2.16 


34 


17 32 54 


-2512 25 




255 


244 


168.9 


3.17 


236 


7.98 


1.38 


169.9 


2.54 


3.91 


1.87 


35 


17 33 01 


-25 46 00 


-89 


133 


130 


171.7 


2.37 


125 


10.83 


1.80 


171.3 


1.95 


4.48 


2.23 


36 


173011 


-26 48 42 




144 


142 


160.4 


3.76 


138 


9.79 


1.62 


160.1 


3.65 


3.24 


1.51 


37 


173118 


-26 29 36 


66 


111 


111 


169.8 


3.40 


105 


13.92 


2.29 


170.6 


3.11 


4.53 


2.15 


38 


17 32 27 


-2615 49 


74 


127 


127 


172.6 


3.36 


125 


15.51 


2.85 


173.4 


3.26 


4.48 


2.11 


39 


17 38 56 


-24 08 57 


151 


249 


245 


6.5 


6.07 


233 


3.87 


1.04 


7.8 


5.83 
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Identification from fAlves et al7| ( |2QQ7| ) 
Obtained by a Gaussian fitting to the distribution of observed polarization angles (measured from the North Celestial Pole to East). 
Number of stars passing the selection criteria which were used to estimate the mean values, (P) and #<p), (see text). 



similar to the one used by |Pereyra fc M agalhaes| ( |2QQ7[ ) , 
that is, to improve the precision ol the mean values, we 
selected only those objects with observed polarization 
angle 6 obs within the interval (<9 mea n - ^std < obs < 
Omean + ZcFstd) where, d mean and cr std are the mean polar- 
ization angle and standard deviation of each field sample 
(columns 7 and 8, respectively, in Table|3|. The mean 
Stokes parameters, (Q) and (t/), for each neld, were es- 
timated from the individual values for each star 
weighted by the error (c^) according to 



The estimated mean polarization value (P) and its asso- 
ciated error SP are then given by 

(P) = V(Q) 2 + (^) 2 , 

Q\{Q)\ + 5U\(U)\ 



5P 



(P) 



(Q) 



(U} = 



EfaM 2 ) 

E-r 2 
E-r 2 



where 8Q and 8U are the estimated standard deviation 
for the mean Stokes parameters (Q) and (U), respec- 
tively. The mean polarization position angle 0^ is given 
by 

V)=0.5tan- (|T 



The number of stars passing the 2a st( i filter, the ob- 
tained mean polarization, its estimated uncertainty, and 
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Fig. 4. — Representation of the mean polarization degree (filled circles) and dispersion of polarization angles (open circles) for the observed 
areas. The size of the symbols are proportional to the scale indicated over the left-hand corner of the image. The anti-correlation between 
dispersion in polarization angle and mean polarization is clearly seen. 



the polarization angle for the mean polarization vector 
are given in columns 9, 10, 11, and 12 of Table[3j re- 
spectively. Column 13 shows the dispersion of polar- 
ization angles corrected in quadrature by the mean er- 
ror of the polarization position angle, that is, AO = 
{o~ 2 std — (crfl) 2 ) 1 / 2 , where the mean error, (00), was esti- 
mated from {(Jq) = ^o-Qi/N, where gqi is the estimated 
uncertainty of the star's polarization angl^] 

The global polarimetric propertie s of the Pipe nebula 
were already presented in Paper I[ and show some in- 
teresting results. For instance, the obtained mean po- 
larizations for the region of B 59 and a long the stem are 
typical for star for mation regions (e.g., |VrbaeTaL] [ 1993; 



Whittet et aT1|1994 



2001|, while the values obtained lor 
the bowl are unusually high. Another noteworthy result 
presented in Paper I is the apparent general tendency 



of decreasing dispersion in polarization angles along the 
filamentary structure of the Pipe nebula from B 59 to- 
ward the bowl, while the mean degree of polarization in- 
creases toward this end. This effect is better visualized 
by inspection of the image presented in Fig.|4j where we 
represented the obtained mean polarization degree and 

2 The uncertainty of the polarization angle is estimated by error 
propagation in the expression of the position angle 0, which yields 
(jq = \ <jp/P, in radians, or gq = 28?65<7p/P (see for instance, 
|Serko wski 1974) when expressed in degrees. 



dispersion of polarization angles by filled and open cir- 
cles, respectively, scaled to the values of these observa- 
tional quantities. In fact, this figure is more in structive 
than the diagram introduced in Fig. 2 of Paper I[ because 
in addiction to the above mentioned anti-correlation be- 
tween polarization degree and dispersion of polarization 
angles seen along the main axis of the complex, one can 
also see how these two quantities distributes spatially 
over the cloud. In general, one see that fields toward 
lower infrared absorption have the tendency of presenting 
larger values of dispersion of polarization angles. A no- 
ticeable exception is Field 26, located close to the center 
of the area displayed on Fig.[4j which presents the second 
largest dispersion value in our sample (AO = 32? 5). 

All fields, but three, having a rather broad distribution 
of polarization angles (AO > 10°) are in the vicinity of 
B 59. The exceptions are: Field 15, laying in the stem 
almost middle way from B59 to the bowl, the already 
mentioned Field 26, and Field 27, both located in the 
eastern side of the stem, close to the bowl. 

3.2. Polarization maps 

It is instructive to analyze the obtained polarization 
for each CCD field. In Figs.[5]to|9]the obtained polariza- 
tion is overlapped onto the dust extinction maps of the 5 
large areas demarcated in Fig. [I] which cover all observed 
CCD fields except Fields 39 and 45. The histograms give 
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Fig. 5. — Polarization map for Fields 01 to 14 overlapping the dust extinction map of the corresponding area ( |Lombardi et al.|2006| >. The 
overplotted contours are for Ay =2,4 and 8 mag. The length of the vectors correlates linearly with the degree ot polarization according 
to the scale indicated over the left -hand corner of the image. The vertical dashed-line demarcates the limits between the stem (left) and 
B 59 (right), as denned in |Paper~7| Histograms for the distribution of the polarization angles are shown individually for each field. The 
identification of the fields is given in the upper right-hand corner of the histograms. The o verplotted gaussian c urves are for comparison 
purposes only. The 'star' symbols indica te the location of the identified candidate YSOs by |Forbrich et alT] |2009). Note that the location 
of the young stellar cluster identified by |Brooke et al.| ( |2007| in the heart of B 59 was omitted. The source on the west side of B 59 core 
is the Herbig Ae/Be star KK O ph which is very likely associated to the cloud. The two other objects at the east of B 59 are sources 11 
(north) and 16 (south) listed by |Forbrich et al.| ( |2009) . 



the distribution of obtained position angles for each field, 
identified in the upper right corner. The gaussians rep- 
resent the best fit to the distribution and are showed for 
comparison purposes only - they help us to visualize how 
the distributions of some fields d epart from the "normal 
distri bution. In a classical work, [Chandrasekhar & Fermi 



line-of-sight toward one of the m ost opaque regions of 
the entire n ebula, the B 59 region ([Roman- Zuhiga et al.| 
2007 2009), is the observed field with the smallest num- 



(1953) obtained a reasonably accurate estimate for the 



field strength in the diffuse ISM by directly relating the 
dispersion in polarization position angle to the ratio of 
two energy densities: the energy density of the uniform 
component of the field and the energy density of turbu- 
lence. Since then, it is widely accepted that the mean 
value of the distribution of polarization position angle 
obtained from a polarization map gives the angle of the 
mean or uniform (large-scale structured) magnetic field 
for the region under investigation, while the dispersion in 
the distribution gives information about the statistically 
independent non-uniform (turbulent or random) compo- 
nent of the magnetic field (a de tailed discussion concern- 
ing this subject can be found in M yers fc Goodman) 1991 , 
and references therein). 

The effects of the high interstellar absorption in some 
of the observed fields are clearly seen on the distribution 
of the measured stars. For instance, our Field 03, with 



ber of stars with P/<Jp > 5 (21 stars only). Its histogram 
of observed polarization position angles and the obtained 
mean position angle, 6/ P ) = 48? 8 (Table|3|, indicate that 
most of those stars belong to the right-nand tail of the 
polarization position angle distribution given in Fig. [3] 
(top panel). Although the obtained large dispersion or 
position angles - which is due to 6 stars - the distribu- 
tion of the remaining stars is rather narrow, as seen in 
the histogram for Field 03 shown in Fig. [5] The two po- 
larization angles on the right-hand side of the main dis- 
tribution (0 = 64?9 and 70? 1) correspond to [BHB2007] 2 
and [BHB200 71 1, respectively, supposed to be c andidate 
young stars ( jBrooke et al.||2007| [Fbrbrich et al.||2009[ ). It 
is noteworthy that high resolution optical images of the 
region show a "light cone-shaped" which apparently em- 
anates from these stars and illuminates the surrounding 
dust material. Interestingly, both observed polarization 
vectors are almost perpendicular to the symmetry axis 
of this cone. 

From the histograms shown in Fig. [5] we note that most 
of the fields presenting large dispersion of polarization 
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angles suggest a multicomponent structure, in special, 
Field 06 presents a very interesting geometry for the ob- 
tained distribution of the polarization vectors and de- 
serves further comments (see §[5|. Fields 01 to 04 show 
distributions of polarization angles with many stars hav- 
ing values between 0° and 40°, while the remaining fields 
given in Fig.[5j already show distributions with polariza- 
tion angles between 160° and 20°, likely what was ob- 
tained for most of the other observed fields. 

The area covered by Figj5] contai ns, apart from the 
young stellar cluster identified by Brooke et al. (2007) 
embedded in B 59, three of the six candidate YSOs found 



by Forbrich et al. (2009). One of them is the well known 
KK Oph, a pre-mam-sequence binary with 1.6" separa- 
tion and suggested to constitute a H erbig Ae star with 
a classical T Tau ri companion (e.g., |Herbig||2005| [CaF 



mona et al. 2007, and references therein) . Althougn cle 



Winter &; The| (|1990) attribute a distance of 310 pc to 
this star, it is commonly ac cepted a distance of 160 pc 



(e.g., Hillenbrand et al.|1992 ) suggesting that this object 
may have been fo rmed by materia l form erly associated to 



the Pipe nebula. Car mon a et al.| (120071) estimate an age 
of about 7Myr to this system, tnat is, from 5 to 6 Myr 
older than the estimated age of the YSOs in B 59. The 



two other objects are sources 11 and 16 in the Forbrich 



et al. (2009) candidate YSOs list, and have lines-of-sight 



toward Fields 09 and 11, respectively, close to the tran- 
sit ion betw een the B 59 and the stem regions, as defined 
in P aper I| These sou rces were spectroscopically studied 
by Covey et al. (2010), who confirmed the youthful char- 
act ieTTTtEeTatter, and found that it is a visual binary, 
while the former presents an ambiguous spectra, that is, 
it may either be a rather young object or a reddened 
giant /subgiant. 

Figure[6] displays the middle portion of the stem. We 
observed seven fields in this area. The obtained his- 
tograms seem to present a kind of transition between 
the characteristics observed for the B 59 region and the 
ones for the bowl. That is, Field 15 (one of the fields with 
AO > 10°) have a distribution that resembles the ones 
obtained for the fields in B 59, while very close to it one 
see Field 16 which shows a distribution with a dispersion 
typical of the ones presented by fields in the bowl, how- 
ever, centered around 0°. On the left-hand side of this 
figure there is Field 22 showing polarization properties 
with all the characteristics observed in the bowl, that is, 
low dispersion of polarization angles and centered around 
160°. Two of the |Forbrich et al.j ( [2009] ) YSO candidates 
are located in this area. One of them very close to the 
center of our Field 17 (source 26), the other one close to 
the border of Fie ld 18 (source 24). N one of these sources 
were studied by Covey et al.| ( |2010[ ), who on the other 
hand, investigated two other sources that turned out to 
be OH/IR stars, likely residing in the Galactic Bulge. 

Figure[7| displays the "transition region" between the 
stem and the bowl, as denoted by the vertical dashed-line. 
In this area we find two of the three fields presenting 
broad distribution of polarization angles not belonging 
to the B 59 vicinity, Fields 26 and 27, the former shows 
a distribution of polarization vectors that resembles the 
one observed for Field 06, suggesting that there may be 
some similarities between the physical properties of both 
cores, while the distribution for the latter clearly shows 
a bimodal distribution of polarization vectors. All fields 



presenting particularly interesting polarization distribu- 
tion are separately discussed in §[5] The four eastern 
fields of this area present polarimetric characteristics of 
the bowl, that is low dispersion of polarization angles, 
rather high polarization degree and a mean polarization 
angle centered around 160°. A detail that calls our at- 
tention is the polarization probed by Field 20 (upper 
right corner of the figure). While all other fields shown 
in this figure present a distribution of polarization angles 
centered around ~160°, the distribution of polarization 
angles for Field 20 is centered around ~30°. This is the 
second less absorbed field (Ay ~ 1^5), so that the po- 
larization mapped by these stars may be mainly caused 
by a background medium. 

The area displayed in Fig. [8] is located to the north of 
the one displayed in Fig. [7] and covers mostly the more 
diffuse medium of the Pipe nebula, except for Field 29 
with line-of-sight toward a portion with higher extinc- 
tion. Although this field presents a rather large disper- 
sion of position angles, as compared to the other fields in 
the bowl, its rather high mean polarization and mean po- 
larization angle centered around 160°, are characteristics 
of that part of the complex. 

Figure|9] displays eleven fields observed in the bowl 
area. The main characteristics of the fields observed to- 
ward this region of the Pipe nebula are the high degree of 
polarization and the highly aligned polarization vectors, 
as testified by the low dispersion of polarization angles 
shown by the histograms displayed on this Figure. 

It is also interesting to analyse the distribution of mean 
polarization angle as a function of the right ascension of 



the observed fields. Such distribution is shown in Fig. 10 



and as already mentioned, due to a fortuitous coinci- 
dence this celestial coordinate correlates quite well with 
the field's position along the long axis of the Pipe nebula. 
Most of the obtained mean polarization angles are in the 
interval 0^ ~ 180° ± 20°, indicating that the local uni- 
form magnetic field is somewhat aligned perpendicularly 
to the main axis of the cloud complex. 

Apart from four fields, identified in Fig.[l0j the distri- 
bution of the remaining mean polarization position an- 
gles seems to follow a pattern. The values obtained for 
fields toward directions having lower infrared absorption, 
represented by open dots, present a rather constant value 
all over the stem of the Pipe, including the B 59 region, 
except for two of the four mentioned fields (Fields 20 
and 27). On the other hand, the distribution shown by 
fields with rather large infrared absorption, represented 
by filled dots, is more interesting. Again, apart from 
the other two identified fields (Fields 03 and 06, both 
in the B 59 region), which as we have mentioned earlier 
show some kind of peculiar characteristic, one see that 
the mean polarization angles for these fields seems to be 
rather constant (0(p\ ~ 180°) from B59 to almost the 
center of the stem, then decrease slowly until close to our 
arbitrary border of the bowl region, and rise up again by 
a small value and became almost constant (0(p) ~ 170°) 
in the bowl. 

This behaviour somehow suggests that the uniform 
component of the magnetic field is "uniform" in the sur- 
rounding diffuse medium but presents small systematic 
variations along the dense parts of the complex. A re- 
markable point to be noted is the fact that the right as- 
cension of Fields 20 and 27 somehow coincides with the 
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one where the mean position angle of the fields associated 
to cores seems to reach its smallest value. Unfortunately 
our observational data do not allow us to investigate fur- 
ther this coincidence. 

3.3. Deriving Ay from 2MASS data 

It is instructive to compare the obtained mean polar- 
ization with the interstellar extinction acting on each ob- 
served line-of-sight. The mean extinction in each field 
could have been e stimated by the use of the extinction 
map obtained by Lombardi et al. (2006). In fact, we 
started using their image with this purpose, however, 
as already mentioned, some of the observed fields con- 
tain areas of high interstellar absorption that were not 
probed by our stellar sample. Thus, simply averaging the 
infrared extinction over the observed field would provide 
a larger value for the reddening than the one actually 
probed by our obsereved stars. Because of that, we de- 
cided to use another approach, that is, the mean extinc- 
tion in each field may be determined by assuming that 
the old bulge population present in each observed vol- 
ume has an upper giant branch similar to that found in 
K, J — K color magnitude diagram s (CMD) of Baade's 



diag: 

window (see e.g.,|Frogel et al."|1999|). We proceeded as- 



suming that the upper giant branch in each of our ob- 
served fields is comparable to and has the s ame slope as 
the ex tinction-corrected template derived by |Dutra et al. 
(2002), given by 



(K s )o = "7.81 x (J - K s )o + 17.83 



(1) 



This assumption is perfectly justified, because those au- 
thors applied this template to study the interstellar red- 
dening in a volume that partially contains the Pipe neb- 
ula. We also assumed that the relation between extinc- 
tion and reddening is given by 



A Ks = 0.670E(J - K s ). 



(2) 



From the (Ks, J — Ks) CMD values of each star in the 
field, we calculated the shift along the reddening vector 
given by Eq. ^ required to make it fall onto the ref- 
erence upper giant branch, Eq. (fil. Since the adopted 
template appropriately describes the upper giant branch 
locus for stars with 8 < (Ks)o < 12.5, all star present- 
ing a corrected Ks magnitude outside this range were 
excluded from our mean absorption estimate, and sim- 
ilarly to what we have done when estimating the mean 
polarization, a 2-a filter was applied to the obtained dis- 
tribution of E(J — Ks) and the field's extinction value 
was taken as the median of the distribution of stars pass- 
ing the clipping selection. The estimated mean Ak s val- 
ue s were then convert ed to Ay by the relation derived 
by |Dutra etal] (|2002|) , i.e. A Ks /A v = 0.118. To illus- 
trate the method used to estimate the mean interstellar 
absorption, we show, in Fig. [II] the CMD obtained for 
our Field 43. It is clearly noticeable that most stars 
brighter than Ks ~ 12 m , in this field, are reddened by 
about E(J — Ks) = 0™5. Stars used in our estimate 
of the mean interstellar absorption are represented, in 
Fig. [IT] by filled circles. In general, an analysis of the 
(J — H, H — Ks) color indices shows that most of the 
stars fainter than Ks ~ 12 m in our diagrams are likely to 
be main-sequence stars of earlier types (typically B-G). 



The top-panel of Fig. 12 shows the obtained CMD for 
all observed star, in our sample, with P/ap > 5 and 
identified in the 2MASS catalogue. For comparison, the 
(J — H) — (H — Ks) diagram for the same stars is given 
in the bottom-panel. As one can see, most of the ob- 
served stars seems to occupy the area corresponding to 
normal reddened stars. The stars in this zone could also 
be dereddened onto intrinsic color lines by extrapolation 
from the observed stellar locus along the appropriate red- 
dening vector. 

The obtained values of Ay and their estimated uncer- 
tainties are given respectively in columns 14 and 15 of 
Tabled The later were estimated from the standard de- 
viation of the obtained distribution of E(J — Ks), before 
applying the 2-a filter, being that the stellar photometric 
errors have not been taken into account. 

Although the extinction can reach very high values 
toward some of the observed cores, one see in Table[3] 
that, as expected, our optical polarimetric survey is prob- 
ing the less absorbed areas only (e.g., from Ay > 0™6 
to Ay < 4^6) . It is impor tant to note that although 
|Roman-Zuhiga et al.| ( |2007[ ) found evidences that the 
extinction law prevailing m the densest regions of B 59 
agrees more closely with a dust extinction with a total 
to selective absorption Ry = 5.5 we adopted the typical 
values for the diffuse interstellar medium, since we are 
studying regions with extinctions well below Ak < 2 m . 

As an independent control of the method used to eval- 
uate the interstellar absorption towards the observed ar- 
eas, one may compare the results obtained for our ex- 
ample field, Field 43, with the ones we would have ob 
tained from the extinction map produced by 



Lombardi 



20061). As one may verify from visual inspection 
our stellar sample covers rather uniformly the 
area deEn ed by Field 43 w hich means that, simply aver- 
aging the |Lombardi et al.| extinction over this area will 
probably yield values that are representative of the mean 
extinction probed by the observed stars. Such procedure 
give us Ay = 3™28 ± i™26, which agrees quite well with 
the value given in Table[3J 

4. POLARIZING EFFICIENCY TOWARD THE PIPE 
NEBULA 

The degree of polarization produced for a given amount 
of extinction is referred to as the polarization efficiency 
of the intervening dust grains. This efficiency of polar- 
ization depends on both the nature of the grains and 
the efficiency with which they are aligned in the line- 
of-sight. The most efficient polarization medium con- 
ceivable is obtained by modeling the dust grains as in- 
finite cylinders (very long in comparison to their radii) 
with diameters comparable to the wavelength, perfectly 
aligned with their long axes parallel to one another and 
perpendicular to the line of sight. For such a model, Mie 
calculations for particles with dielectric optical proper- 
ties, place a theoretical upper limit on the polarization 
efficiency of the grains due to directional extinction at 
visual wavelengths of p/Ay < 14% mag -1 (see, for in- 
stance, |Whittet 2003). The observations, however, show 
that the upper limit predicted by this very idealized sce- 
nario is far from being reached. In general, studies of in- 
terstellar polarization demonstrate that the efficiency of 
the real Galactic interstellar dust as a polarizing medium 
is more than a factor of 4 less than predicted for the ideal 
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Fig. 9. Same as Fig.[5]for Fields 33 to 35, 37, 38, 40 to 44, and 46. 
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Fig. 10. — Distribution of mean polarization angle, #(p), as a 
function of the right ascension of the observed field, which cor- 
relates quite well with its position along the long axis of the Pipe 
nebula. Filled (blue) and open (red) dots represent values for fields 
associated to lower and higher infrared absorptions, respectively. 
The former are mostly fields outside the main structure of the Pipe 
nebula, namely, Fields 02, 04, 07, 08, 16, 19 to 21, 23 to 25, and 
27. The gray bars give the interval defined by 0^ ± AO, where 
AO is the dispersion of polarization angles (see Table[3}. 



polarizing medium. The observational upper limit on the 
ratio of polarization to extinction for the diffuse interstel- 
lar medium is given by p/Ay ~ 3% mag -1 ( |Serkowski| 
eTaI1[T975 ). 




Considering that our sample contains a rather large 
number of objects, in the bowl region, showing outstand- 
ing degrees of polarization, it is natural that we try to 
investigate if these observed areas present unusual polari- 



FlG. 11. — Color-magnitude diagram for stars in Field 43. Stars 
passing the selection criteria (see text) and used to estimate the 
field's mean interstellar absorption are represented by filled (red) 
circles. The straight line represents the reference upper giant 
branch (equation [lj. It is clear the gap between the standard 
reference line and observed stars brighter than Ks ~ 12 m , which 
corresponds to an interstellar reddening of about E(J—Ks) = 0™5. 
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Fig. 12. Top panel: K S -{J - K s ) CMD for all observed star 
with P/crp > 5 identified in the 2MASS catalogue. The straight 
line represen ts the reference upper giant branch, EqJ^ obtained by 
|Dutra et alT] ( |2002| ) . Botton panel: (J — H) — (H — K s ) color-color 
diagram tor trie same stars. The theoretical locus for the main- 
sequence stars is represented by the continuous line, while the white 
(red dashed-line) represents the giant branch stars. The absorption 
vector indicated in both diagrams corresponds to Ay = 5 m . 

metric properties as compared to the common Galactic 
interstellar medium. 

The diagrams shown in Fig. were constructed in or- 
der to investigate the obtained ratio between our esti- 
mated mean degree of polarization and mean total inter- 
stellar absorption for the observed fields toward the Pipe 
nebula — error bars were omitted in these diagrams for 
the sake of clarity. The plot of mean polarization versus 
total visual absorption given in the Fig.[l3] (top panel) 
shows that basically all data points lie on or below the 
line representing the usual relation p/Ay « 3% mag -1 
— the two points appearing above this line represent 
data obtained for Fields 38 and 41, however, taking into 
account the estimated 1-cr uncertainties for the mean 
degree of polarization and total interstellar absorption, 
these two fields may also obey the above relationship. 
On the one hand, this result indicates that the interstellar 
material composing the Pipe nebula follows the usual be- 
haviour of the common diffuse interstellar medium. On 
the other hand, as one can see from the values tabu- 
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Fig. 13. — Top panel: Plot of mean polarization (P) versus to- 
tal visual absorption (Ay) derived from the 2MASS data for the 
observed stars with P/crp > 5. The solid line represents optimum 
alignment efficiency (P(%) = 3 x Ay). Middle panel: Polariza- 
tion efficiency (P/Ay) versus visual absorption Ay. Bottom panel: 
Distribution of the polarization efficiency as a function of the right 
ascension of the observed field. Symbols have the same meaning 
as in Fig.[l0l 

lated in Table[3] we have found levels of mean degree of 
polarization that are unusual for the same interstellar 
material. 

Several previous investigations have suggested that the 
polarizing efficiency of the interstellar dust declines sys- 
tematically with total extinction, as one probes progres- 
sively denser environments within a dark cloud (e.g., 
Goodman eFaLl|1992| [1995] |Gerakines et al.|[l995l ). The 
obtained diagram of polarizing efficiency, p/Ay, as a 
function of the interstellar absorption (Fig. [13] - mid- 
dle panel), does not show clearly this tendency, at least 
not for the covered interval of interstellar absorption. 
In fact, on the contrary, if we exclude Field 02, which 
shows the lowest interstellar absorption and a polariza- 
tion efficiency of almost 3 % mag -1 , the other fields show 
a tendency of increasing efficiency with the interstellar 
absorption. 
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Fig. 14. — Polarization vectors overlaid on the optical image of 
Field 06. All measured polarization for this field are represented in 
this figure and not only the ones having P/crp > 5. The observed 
orientation of the polarization vectors seems to embrace the dust 
core whose existence is suggested by the scarcity of observable stars 
to the left of the center. The length of the vectors correlates linearly 
with the degree of polarization according to the scale indicated over 
the left-hand corner. 



More interestingly is the diagram shown in the bot- 
tom panel of Fig.[l3 , which shows the distribution of the 
estimated polarization efficiency of the observed fields 
as a function of their position along the long axis of 
the Pipe nebula. It is known that variations in polar- 
ization efficiency might result from changes in physi- 
cal conditions that affect alignment efficiency, such as 
temperature, density and magnetic field strength, or in 
grain properties such as their shape and size distribution 
and the presence or absence of surface coatings. Most 
of the observed fields in the stem (including its tip — 
the B 59 region) , present a polarization efficiency around 
p/Ay ~ 1 % mag -1 , then it rises up and down when one 
move along the bowl from west to east, reaching values of 
about p/Ay < 4 % mag -1 . Summarizing, although show- 
ing an interesting behaviour, the global properties of the 
probed dust material composing the Pipe nebula does 
not seem to present any special peculiarity, when com- 
pared to the common diffuse interstellar medium, that 
could explain the observed high degrees of polarization. 
However, one notice a clear difference between the be- 
haviour shown by the polarimetric properties presented 
by fields located in the stem and in the bowl. 

Although the division between the regions denomi- 
nated B 59 and stem was chosen rather arbitrarily (in 
|Paper I[ it is characterized by a rising on the degree of 
polarization), one notice an interesting feature in the bot- 
The polarization efficiency seems 



torn panel of Fig. [13 
to increase along the dust filaments probed by our sam- 
ple when we move from the B 59 region to the stem. It 
happens only for the fields of the stem shown in Fig.[5j 
after that, the ratio p/Ay returns to the typical value of 
~l%mag -1 observed for B 59 and the remaining fields 
in the stem. This behaviour can be an indication that 
distinct physical regimes may be acting on different frag- 
ments of the stem. For instance, variations of the value of 
p/Ay may arise where the local magnetic field is not or- 
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Fig. 15. Same as Fig.|14|for Field 26. 

thogonal to or its direction varies along the line-of-sight, 
or where the processes responsible for grain alignment 
change for some reason (grain composition, size, shape, 
etc). The interested reader will find a good r eview on 
the ef f icienc y of grain alignment in the work by | Whittet] 



et al. (2008, and references therein). It is worthwhile to 
mention that the point where the polarization efficiency 
returns to its typical value of ~l%mag -1 almost coin- 
cides with the place where one noticed the value of the 
mean polarization angles started decrease (see Fig. 10). 

All these results reinforce once more how interesting is 
the Pipe nebula and suggest that this complex may be 
a testbed for different theories of dust grain alignment 
efficiency. 

5. FIELDS SHOWING INTERESTING POLARIZATION 
DISTRIBUTIONS 

Inspection of Figs. [5] to [9] shows that some of the ob- 
served fields present remarkable polarization geometries. 
For many of them one clearly note that the obtained 
polarization angles for the objects in the field suggest 
a multicomponent, or in some cases a hoop-like, distri- 
bution. As one have seen, in most cases the mean po- 
larization vector is aligned perpendicularly to the long 
axis of the Pipe nebula, but there is the case of Field 20 
(see Fig.[l0|) where the distribution of polarization angles 
does not follow the average behaviour for the region. Be- 
low we introduce three of the most interesting observed 
fields, and comment on the fields having high mean po- 
larization ((P) > 10%). 

5.1. Field 06 

The polarization map for Field 06 is shown in Fig. [14} 
one of the most interesting distribution in our survey. 
This is one of the four fields we have identified in Fig. [l0| 
as having a mean polarization angle which seems to dis^ 
agree from the pattern observed for the cloud complex. 
In order to emphasize the geometry of the magnetic field 
in this region, all star for which polarization has been 
measured are represented in the map. The polarization 
vectors seem to suggest that the local magnetic field fol- 
lows the border of the dust cloud evidenced by the higher 
interstellar absorption noticed to the left of the center. 
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Has this core been modeled by the field or, on the con- 
trary, was the field shaped by the core? In any case, this 
seems to be an interesting region which deserves further 
investigation. 

5.2. Field 26 



The polarization map obtained for Field 26, Fig. 15 
seems to be the result of a mixture of two distributions, 
a main component centered around 160° (see also the his- 
togram introduced in Fig. [7]), combined with an hoop-like 
component. An interesting point is that, as one may ob- 
serve in the polarization map, the surveyed area seems to 
show different characteristics toward directions located 
northern and southern of the densest parts of the cloud 
— visually characterized by the absence of stars. Ap- 
parently, at the south only the main component of the 
distribution (6 ~ 160°) is present, while at the north we 
observe the presence of both distributions. An inspection 
of Fig. [7] shows that the northern part of this field probes 
a more diffuse part of the interstellar material, as what 
happens in the case of Field 27 (see below), while the 
southern stars have line-of-sight toward a volume pre- 
senting higher ex tinction. One of the cores studied by 



Frau et al. 



_ ( |2Q1Q[ ), who used the IRAM 30-m telescope 
to carry out a continuum and molecular s urvey toward 
four o f the starless cores from the list of lAlves et al.l 
(2007), is Core 48, which is associated to the higher in- 



terstellar absorption shown in Fig. 15 The radio data 



indicates that, although being very diffuse, this core has 
a strong dust emission, and their molecular analysis sug- 
gests that chemically it seems to be in a very early stage 
of evolution. 

5.3. Field 27 

There is no dense core associated to the volume probed 
by this field, and it is other of the fields having mean po- 
larization angle not fitting in the main pattern of mean 
position angles, as defined in Fig. fTol The distribution of 
polarization vectors shown in Fig. |16|( see also the histo- 
gram of polarization angles shown in Fig. [7]) clearly shows 
a bimodal distribution with mean angles values centered 
on ~135° and ~155°. Both components seem to be well 
distributed all over the surveyed field. 

5.4. Distribution of polarization and position angles as 
function of the 2MASS Ks magnitude for Fields 
26 and 27 

The top panels of Fig. [IT] display the measured polariza- 
tion angles as function of the 2MASS Ks magnitude for 
Fields 26 and 27. An interesting result comes out from 
these diagrams. One clearly notice that the distribution 
shown by Field 27 (right panel) is rather defined by the 
stellar Ks magnitude and occupies different regions of 
the diagram. Stars having Ks > 12 m , that is statisti- 
cally popu lated by main sequence stars, as already men- 
tioned in § |3.3[ are mainly associated to the component 
having higher mean angle, while stars having Ks < 12 m , 
statistically populated by giant stars, are basically asso- 
ciated to the component having lower mean angle. This 
is indicated by the horizontal and vertical dashed lines 
positioned at 9 = 150° and Ks = 12^0, respectively. 

The distribution presented by Field 26 is rather differ- 
ent but shows some of the characteristics presented by 
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Fig. 16. Same as FigTuflfor Field 27. 
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Fig. 17. — Distribution of polarization degrees (bottom panels) 
and polarization angles (top panels) as a function of the K magni- 
tude. Data obtained for Field 26 is represented on the left panels 
and for Field 27 on the right panels. All observed stars in each field 
were used to construct these diagrams. The horizontal and vertical 
dashed lines, represented in both top panels, were arbitrarily posi- 
tioned at = 150° and -f^MASS = 12™0, respectively (see text). 
Mean uncertainties of the quantities are indicated by the horizontal 
and vertical bars on the lower left corner of each diagram. 



Field 27. For the sake of comparison, we have repre- 
sented the same horizontal and vertical dashed lines in 
both diagrams. While the polarization angles observed 
for Field 27 are restricted between - 120° and 170°, 
Field 26 presents basically all values of polarization an- 
gles. However, as observed for Field 27, most of the stars 
in Field 26 fainter than Ks = 12 m has polarization angle 
larger than ^140-150°, suggesting that the same kind of 
interstellar structures may be present toward both line- 
of-sights, which are separated about 20' from each other. 

It is also interesting to compare the distribution of de- 
gree of polarization as a function of the stellar magni- 
tude (bottom panels). First of all, one notices that al- 
though Fig. [7] seems to indicate that the line-of-sight to- 
ward Field 27 is less affected by interstellar absorption 
than Field 26, the measured polarization for the latter 
is generally smaller than the one obtained for stars in 
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Fig. 18. Same as Fig.[l4]for Field 35. 

the former field — it must be noted, however , tha t the 
estimated average interstellar absorption in § [3.3| is es- 
sentially the same for both fields (see Tablep]). The 
Ks — (J — Ks) CMD for the observed stars m Field 
27 suggests that the interstellar absorption toward this 
line-of-sight is rather more uniform than the one probed 
by stars in Field 26, as one should expect from the dust 
extinction map obtained by Lombardi et al. (2006) and 
shown in detail by our Fig. [7] Thus, the estimated av- 
erage interstellar absorption ior Field 27 is more repre- 
sentative of what we have all over the surveyed volume, 
while the one estimated for Field 26 is a mean between 
regions showing rather high absorptions, e.g. toward the 
southern area of the CCD field, with regions not so ab- 
sorbed probed by the stars located in the northern area 
of the CCD field. 

5.5. Comments on the Fields with high mean 
polarization degree 

Five of the observed fields present mean degree of po- 
larization (P) > 10%, they all lay in the bowl and are 
Fields 35, 37, 38, 40, and 41. In Fig.|9| these fields are al- 
most aligned along the diagonal crossing the image from 
the upper left-hand to the lower right-hand corner. The 
main characteristics of these fields, apart from the high 
value of observed polarization degree, is the very low dis- 
persion of polarization angles, which suggests that the 
turbulent energy prevailing on the observed cores must 
be quite low (see §|6|. In particular, it is noticeable 
the quite low dispersion presented by Field 35 (see also 
Fig. 18), the lowest in our survey, with a rather "normal" 
distrTBution. 

Although also presenting a very low dispersion, Field 
38, the one with the highest mean polarization in our 
survey, shows a fairly asymmetry in the observed dis- 
tribution of polarization angles. As shown in Fig.[l9j it 
may be caused by two dust cloud components along the 
observed line-of-sight, each one subject to slightly differ- 
ent orientations of ambient magnetic fields. These clouds 
may be associated to the two main velocity components 
that seem to characterize the kinematics of the 'bowl" 
(e.g., |Muench et al.|[2007| ), even though they have not 




170 180 
polarization angle (°) 

Fig. 19. — The distribution of polarization angles for Field 38 is 
clearly asymmetric suggesting two components. Our best fitting 
represented by the full line (red) is the result of two Gaussian 
components, one centered at 170? 2, a st( i = 1?83, and other at 
176? 2, (T st( i = 2? 36, represented by the dashed lines (blue). All 
observed stars in this field have P/crp > 11, being that most of 
them have a much larger signal-to-noise ratio, meaning that the 
theoretical uncertainties of the estimated polarization angles are 
in general much smaller than 2? 6. 

detected two C 18 components toward their observed 
line-of-sight through this field. 

The distribution of polarization angles as a function 
of the 2MASS Ks magnitudes does not present any re- 
markable feature, unless for the fact that the 6 brightest 
stars in the field (Ks < 8™0) have polarization angles be- 
tween 169° and 172? 5, while the remaining stars present 
a rather normal distribution between ~ 165° and 180°. 



Field 40 c ontains other of the cores observed by |Frau 
et al. ( |2010| ), Core 109. The radio data show that this 
object presents a strong dust continuum emission, is the 
densest among the four investigated cores, and one of the 
most massive. The interstellar extinction experienced 
by the observed stars is very nonuniform, ranging from 
A v « 2 m to A v > 5 m . Interestingly the observed 13 CO 
molecular emission shows a double velocity component 
(Alves et al., in pr eparation), which is not seen in C 18 
( Muench et al.|[2007| ratified by the work in preparation 
by Alves et al"), and could explain the asymmetry of 
the distribution of polarization angles which, as observed 
for Field 38, is also noticed for this field but this time 
due to a small excess in the left wingof the distribution 
(see distribution introduced in Fig.[9|. Analyzing the 
distribution of the polarization angles as function of the 
2MASS Ks magnitudes one obtained that this excess 
is due to stars brighter than Ks ~ 11™ 5, which are in 
average more affected by the interstellar absorption and 
present higher mean degree of polarization. Although 
located in the bowl, supposed to be the less evolved region 
of the Pipe nebula, th e molecular investigation conducted 
by Frau et al. ( 2010| ) indicated that the core may be one 
of the chemically most evolved in their molecular survey. 

6. THE STRUCTURE FUNCTION OF THE POLARIZATION 
ANGLES IN THE PIPE NEBULA 

6.1. Basic definitions 

The second-order structure function (hereafter SF) of 
the polarization angles, (AO 2 (I)}, is defined as the av- 
erage of the squared difference between the polarization 
angles measured for all pair of points separated by a dis- 
tance I (e.g. see equation 5 of Falceta-Goncalves et al. 
2008). Thus, the SF give information on the behavior 
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Fig. 20. — Square root of the second order structure function of the polarization angles, (A<l>(/) 2 ) 1 / 2 , for all the individual fields (except 
for Fields 3, 6, and 26) observed in the Pipe nebula. The units of (A<l>(/) 2 ) 1 / 2 are degrees. The upper row panels have a larger range of 
values in the abscissa axis. Fields 1 to 8 are located in B59. Fields 9 to 24 and 27 are in the stem. The rest of the fields belong to the 
bowl. The dashed line show s the 16? 6, which is the value of the dispersion of the position angles for turbulent and magnetic equipartition 
( |Troland fe Crutcher|2008] ). 



of the dispersion of the polarization angles as a function 
of the length scale in molecular clouds. Recently, it has 
been used as a powerful statistical tool to infer infor- 
mation of the relationship between the large-scale and 
the turbulent components of the magnetic field in molec- 



ular clouds (TFalceta-Gon calves et al. 2008 Hildebrand 
~HT" 



et al.||2009[ 

cal sample ol the polarization data in the Pipe nebula, it 
is interesting to compute the SF along the Pipe nebula 
to scales up to few parsecs. For a qualitatively discus- 
sion we will first use the square root of SF, also called 
the a ngular dispersion function or ADF (Poidevin et al. 
20101). The use of the ADF instead of the ^' allows a 



is the highest. A steep slope is an indication that the 
large-scale magnetic field orientation in the plane of the 
sky changes significantly. Fields 3 and 6 have dispersion 
values at scales larger than ~ 0.1 pc close to the expected 
maximum dispersion that would be obtained in case of 
a purely random pol arization angle d istri bution, ~ 52° 



oude et al.||2009p . Given the large statisti- ( Poidevin et al.|2 010). As pointed in § 5.1 for Field 6 this 



more straightforward comparison of the behavior of po- 
sition angle dispersion as a function of the length scale. 
Then, we use the SF to c ompare our statistical sample 



with the previous works (Falceta-Goncalves et al. 2008 
Houde et al]|2009| ). 



6.2. Qualitative analysis 

We have first computed the ADF for all the individual 
fields using a logarithmic scale between 8" (5.6 mpc) and 
11. 7 8 (0.35 pc). This range was selected in order to have 
a good statistical sample. Figure[20| shows the ADF for 
all the Fields but 3, 6, 26. These three fields, which 
are the ones that exhibit the highest polarization angle 



21 



dispersion (see Table pi), are shown separately in Fig. 
There is a clear trena in the distribution of the AlTT 
along the Pipe nebula (see Fig. 20). On one hand, fields 



in B 59 (1-8) not only show a higher polarization angle 
dispersion at all the observed scales but the ADF slope 



is due to a strongly distorted field surrounding a core. On 
the other hand, all the Fields in the bowl (28-46) not only 
have a remarkably small dispersion of the position angles 
( Alves et al.||2008[) but this trend is also observed in the 
ADF at all the observed scales. Indeed and in contrast 
with B 59, the almost flat slope of the ADF in the bowl 
Fields indicate that the projected magnetic field in the 
plane of the sky is very uniform. The ADF behavior 
of the stem Fields is intermediate between that of B 59 
and the bowl. However, the global ADF properties of 
field 26 differ from the general trend found in the stem: 
it shows an unusual high dispersion, ~ 40° at all scales. 
Compared with other Fields with also a high dispersion 
(e.g. 3 and 6) the ADF slope of Field 26 is relatively 
flat. Field 27, the one with a bimodal distribution (see 
§5.3), shows a similar ADF behavior to Field 26 but 
a smaller level: ADF ~ 17° at all scales. Because of 
the peculiarities of these two Fields, we treat them as a 
distinctive region in the Pipe nebula for the SF analysis. 
Here on, we call this region as the " stem-bowl border" . 
We also include Field 20 in this region because its average 
mean direction is quite different from the ones of the rest 
of the stem (see Fig. 10) and it is relatively near to Fields 
26 and 27. 
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Fig. 21. 
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Same as Fig.[20]but for Fields 3, 6 and 26, which show a higher dispersion. 
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Fig. 22. — Second-order structure function of the polarization 
angles, (A<E>(7) 2 ), for the five distinctive regions in the Pipe nebula 
(the list of fields associated with each region is give in Table|4| . 
Right and left panels show the smallest and largest scales for each 
region, respectively. The gray (red) histogram for B 59 shows the 
structure function without Fields 3 and 6. The dashed-lines (blue) 
indicates the best fit of equation [3] for distance up to 0.07 pc. 

Figure|22| shows the SF for the four distinctive regions 
within the Pipe nebula: B 59, the stem (except Fields 
20, 26 and 27), the bowl and the stem-bowl border. By 
combining the different fields of each region, the SF can 
be computed at larger scales, up to few parsecs (see left 
panels of Fig. 22). The general trend described in the 



previous paragraph for the ADF also applies for the SF. 
For example, it is remarkable that the bowl shows very 
low SF values up to scales of few parsecs, with position 
angle dispersion lower than ~ 10°. Yet, B 59 and the 
" stem-bowl border" show an abrupt increase of the SF 
at scales larger than 0.3-0.5 pc. For B 59 this is due to 
the high PA dispersion Fields 3 and 6. Indeed, if these 
two fields are excluded, the resulting SF is smoother at 
these scales. 



6.3. Comparison with Houde et al. (2009) 
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Fig. 23. — Plot of the set of the solutions for the St and a\ 
parameters of the SF. The inner and outer contours show the 
63.3% and 99% confidence regions of the x 2 , respectively. 



The larger statistics obtained by dividing the Pipe in 
four region also allows to better sample the smaller scales. 
The right panels of Fig. [22] show that at scales of few hun- 
dredths of a parsec the SF increases linearly with length 
scale. Yet, at scales of < 0.01 pc, the SF drops fast when 
approaching to the zero length scale. This is a clear in- 
dication that we are starting to resolve the correlation 
length scale for the turbulent magnetic field component, 
which is the scale at which the turbulent energy is dissi- 
pated. In order to approximately estimate the turbulent 
length scale we follow the recipe given in the detailed 
analysis carried out by |Houde et al. ( 2009 ) , where they 
assumed a Gaussian form lor the autocorrelatio n func- 
tion for the turbulence. We use equation 44 from Houde| 
et al. ( |2009[ ) taking into account that the effective an- 
gular resolution of the optical polarization data can be 
considered zero. Therefore, in that equation the angular 



resolution te rm is W = 0. Thus, equation 44 from |Houde 
et al. ( |2009| ) can be rewritten as: 
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TABLE 4 

Structure function parameters for the Pipe nebula 



Region 


ao 

/ radian 2 \ 
^ pc ' 


ai 

(radian 2 ) 


St 
(mpc) 


{SB* IB*)* 


Fields 


B59 


0.44 


0.025 


2.1 


0.4 


1-8 


stem 


0.08 


0.021 


4.8 


0.2 


9-19,21-25 


stem-bowl 


0.38 


0.054 


< 2.1 


0.8 


20,26,27 


bowl 


0.01 


0.008 


4.4 


0.1 


28-46 



Estimated for N = 30 (see text). 



(A$(0 2 ) ~ a l + cn[l - e 



-i 2 /2Sh 



(3) 



The first term, ao Z, gives the large-scale magnetic con- 
tribution to the SF (note that we have adopted a linear 
dependence instead of the original I 2 dependence of the 
aforementioned Eq. 44) . The second term corresponds to 
the turbulent contribution to the SF. S t is the turbulent 
length scale and a\ is a function of the large-scale mag- 
netic field strength, £? , the turbulent component of the 
magnetic field, SB t: and of N, the number of turbulent 
correlation lengths along the line of sightj^] 



ai = (2/N) (5B?/B 2 ) 



(4) 



N can also be understood as the nu mber of independen t 
turbulent cells along the line of sight (|Houde et al.|2 QQ9 ) . 
We have used equation[3] to fit the SF data in the tour 
Pipe regions for the scale range shown in the right panels 
of Fig. [22] These are the scales in which the large-scale 
magnetic field contribution to the SF is basically linear. 
For each of the four Pipe regions a y 2 analysis was carried 
out to find the best set of solutions for the free param- 
eters ao, ai and St. The best-fit solutions obtained are 
given in Table 4] and the 99% and 67% confidence inter- 
vals for a\ ancTot are shown in Fig. [23] The dashed blue 
line in the right panels of Fig. 22 show the best solution 
for each region. We find that the turbulent correlation 
length, St, is in all cases of few milliparsecs. Given that 
the assumption of the Gaussian form for the turbulence 
autocorrelation function is not correct, the found values 
should be taken as an approximation. In any case, the 
right panels of Fig. [22] show that the turbulent corre- 
lation length should be < 0.01 pc. Indeed, at the 99% 
confidence level the x 2 analysis provides an upper limit 
for St of ~ 12 mpc (see Fig. 23). This upper limit is 
slightly lower that the St found lor OMC-1, 16 mpc, from 
submm polarization observations ( Houde et al.| |2009). 

From Eq.[4] we can estimate the turbulent to mag- 
netic energy ratio, SB 2 /Bq, if the number of indepen- 
dent turbulent cells, TV, is known. For optical polariza- 
tion data, the nu mber of independent turbulent cells is 
N ~ A/(V2nS t ) jHoude et al1|2009| ), where A is the 
cloud thickness. A ~ VV(H 2 )/n(H 2 ), where iV(H 2 ) and 
n(H 2 ) are the column and volume densities of the molec- 
ular gas traced by the optical polarization data. 7V(H 2 ) 
can be obtained from the typical visual extinction of the 
observed fields. The average visual extinction for the 
bowl is 3.6 mag, whereas for the rest of the regions is 
2.1 mag. Using the standard conversion to column den- 



3 Note that this term is equivalent to number of the magnetic 
field correlation length along the line of sight introduced by Myers 
|& Goodmanl dl991fc 



sity ( Wagenblast fc Hartquist]|1989| ), these values yields 



to 7V(H 2 ) ~ 4.5 x 10 Zi and 2.6 x 10 21 cm" 2 for the bowl 
and for the re st of the regions, respectively. For the ob- 
served fields, |Paper l| estimated that the volume den- 
sity of the gas associated with the optical polarization is 
n(H 2 ) ~ 3 x 10 3 cm -3 . Therefore, the cloud thickness 
of 0.5 pc for the bowl and of 0.3 pc for the rest of the 
regions. With these values and using for St the range at 
the 67% confidence interval (see Fig. 23 ) we obtain that 
N ranges between 25 and 100 for the Hiowl and 15 to 60 
for the rest of the regions. Nevertheless, a high value of 
N will also reduce significantly the observed polarization 
level. But all the bowl fields and many of the stem fields 
have polarization levels of 4-15% and 3-4%, respectively. 
Therefore, it is unlikely the case of a high TV, at least , 
for these two regions. Indeed, |Myers fc Goodman ( 1991 ) 
estimated that for optical polarizati on observations TV is 
expected to not be larger than ~ 14. |Houde et al. (2009) 
found N ~ 21 for OMC-1 from submm dust polarization 
observations that trace significantly larger column den- 
sities. Therefore, we tentatively adopt a relatively high 
value of N = 30. For this case, the magnetic field appears 
to be energetically dominant with respect to turbulence 
in the Pipe nebula except for the " stem-bowl border", 
where magnetic and turbulence energy appear to be in 
equipartition (see Table [2]). 

(200$ ) 



6.4. Comparison with Falceta-Gongalves et al 



|Falceta-Goncalves et al.| ( |2008[ ) carried out simulation 
of turbulent and magnetized molecular clouds computing 
the effect on the dust polarization vectors in the plane- 
of-the-sky for cases with super-Alfvenic and sub-Alfvenic 
turbulence (i.e., clouds energetically dominated by tur- 
bulence and magnetic fields, respectively). They com- 
puted the SF derived from dust polarized emission as 
well as from optical polarization using background stars 
for the different sub and super-Alfvenic cases, and for 
different angles of the magnetic field with respect to the 
line of sight (see Figs. 6 and 11 of this paper). The SF 
for super-Alfvenic turbulence is clearly higher than the 
one for sub-Alfvenic turbulence: The SF ranges from 0.4 
at the smallest scales up to ^ 1.0 to the highest scales 
(see central panel of Fig. 6 from Falceta-Goncalves et 
al. 2008). For the case of sub-Alfvenic turbulence such a 
high values of the SF are reached only in the cases where 
the magnetic field direction is close to the line of sight. 
For the other cases of sub-Alfvenic turbulence, SF < 0.5. 
Comparin g th e SF obtained in the four Pipe nebula re- 
gions (Fig. 22 ) with the results of Falceta-Goncalves et al. 
(2008) it isclear that B59, the stem, and the bowl are 



compatible with the presence of sub-Alfvenic turbulence. 
The behavior of the SF for the stem-bowl border (SF 
from 0.1 at the smallest scale to > 1.0 at the larger scales) 
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may indicate the case of sub-Alfvenic turbulence with a 
magnetic field near the line of sight rather than super- 
Alfvenic turbulence. Indeed, the only individual field in 
the Pipe nebula that at all scales have a SF compatible 
with the super- Alfvenic turbulence is Field 26. 

6.5. Summary of the SF analysis 

The comparison of the SF derived from our optical 
polarization data with th e ones derived in the works b y 
Houde et al] ( |2009D and |Falceta-Goncalves et al | pQ08 ) ? 
indicated that the Pipe nebula is a magnetically domi- 
nated molecular cloud complex and that the turbulence 
appears to be sub-Alfvenic. Only the region we call the 
stem-bowl border, in particular Field 26, appears to have 
a behavior that is compatible with super- Alfvenic turbu- 
lence. A similar situation seems to apply to the well 
investigated low mass star forming region in the Tau- 
rus complex where there is evidence for a molecular gas 
substrate wi th sub-Alfvenic turbulence and magnetically 
subcritical ([Heyer et"al||2008| iNakamura fc Li||2008| ). 



Hily-Blant & Falgarone (2007) also found that in Taurus, 
the magnetic fields are dynamical important, although 
they found that they are trans- Alfvenic. In addition, 
analyzing the polarization angles at different scales us- 
ing optical and submm obs ervations in several molecular 
cloud yield Li et al. ( 2009[ ) to suggest that these clouds 
are also sub- Alfvenic. 

7. SUMMARY 

The Pipe nebula has proved to be an interesting in- 
terstellar complex where to investigate the physical pro- 
cesses that forestall the stellar formation phases. The 
polarimetric survey analyzed in this work covers a small 
fraction only of the entire Pipe nebula complex, and there 
is no doubt that new data is highly desired in order to 
verify so me of th e speculations settled in this investiga- 
tion. In Paper I, we suggested that the Pipe nebula, a 
conglomerate of filamentary clouds and dense cores, is 
possibly experiencing different stages of evolution. From 
the point of view of the global polarimetric data alone, we 
proposed three evolutionary phases from B 59, the most 
evolved region, to the bowl, the youngest one, however, 
the real scenario seems to be much more comp licated 
than that. As demonstrated by Frau et al. (2010), from 
the point of view of the chemical properties derived for 
four studied starless cores, there does not seem to be a 
clear correlation between the chemical evolutionary stage 
of the cores and their position in the cloud. 

In addiction, the polarimetric analysis conducted here 
suggests that, 

1. Although the unusually high degree of polarization, 
observed for numerous stars in our sample, the 
probed interstellar dust does not seem to present 



any peculiarity as compared to the common diffuse 
interstellar medium. In fact, the fields where the 
high polarization were observed show a polarization 
efficiency of the order of p/Ay « 3 % mag -1 , which 
is the typical maximum value universally observed 
for the diffuse interstellar medium. 

2. Basically all observed fields in B 59 and the Pipe's 
stem present an estimated polarization efficiency 
of the order of p/Ay « l%mag -1 , and all so far 
known candidate YSOs presumed associated to the 
Pipe nebula were found in those regions. 

3. While the value of the mean polarization angle ob- 
tained for fields toward volumes not associated to 
the densest parts of the main body of the Pipe neb- 
ula seems to remain almost constant, the same does 
not happens for fields presenting large interstellar 
absorption, suggesting that the uniform component 
of the magnetic field permeating the densest fila- 
ments of the Pipe nebula shows systematic varia- 
tions along the main axis of the dark cloud com- 
plex. 

4. Analysis of the second-order structure function of 
the polarization angles suggests that in the Pipe 
nebula the large scale magnetic field dominates en- 
ergetically with respect to the turbulence, i.e. the 
turbulence is sub-Alfvenic. Only in a localized re- 
gion between the bowl and the stem turbulence ap- 
pear to be dynamically more important. 
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